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We analyze the propagation of electromagnetic waves in media where the dielectric constants 
undergo rapid temporal periodic modulation. Both spatially homogeneous and periodic media are 
studied. Fast periodic temporal modulation of the dielectric constant of a homogeneous medium 
leads to existence of photonic band-gap like phenomena. In the presence of both spatial and tem- 
poral periodicity the electromagnetic spectrum is described in a four-dimensional cube, defining an 
effective Brillouin zone. In the case of incommensurability between space and time periodicities, 
completely dispersed point spectra exist. 



The advent of materials whose electric permittivity e is 
periodically modulated at the nanometer length scale has 
introduced the important concept of photonic band-gap 
structures. This has in turn made traditional solid-state 
concepts like reciprocal space, Brillioun zones, dispersion 
relations, Bloch wave functions, etc. directly applicable 
to the field of classical electromagnetic wave propagation. 
During the last decades the attractive possible applica- 
tions of such photonic band gap materials have driven 
intense theoretical and experimental studies of the prop- 
agation of electromagnetic waves in^atiallyperiodic and 
disordered dielectric structures 0, This 
has provided the field of photonics with a wide range of 
new applications, mostly related to guided light modes 
0. Recently a magnetic photonic crystal was made 
by periodically modulating the magnetic permeability/i 
of specially constructed material. 

Although there have also been tentative studies of ef- 
fects arising from low frequency 0, and specific cases 
[lH of temporal modulation of the dielectric constant, this 
possibility has been largely ignored. However, as we show 
here, essentially all fruitful concepts from the now ma- 
ture field of photonic band-gap materials can be applied 
in the case of fast temporal modulation of a material's 
dielectric response. Further, the combination of spatial 
and temporal modulation of the dielectric response in- 
troduces intriguing new concepts to the field of classic 
electromagnetic wave propagation. The results are also 
valid for periodical modulation of the magnetic perme- 
ability ji of spatially homogeneous or periodic magnetic 
media. 

First, it is straightforward to realize that temporal 
modulation leads to photonic band structures similar to 
spatial modulation: In a material with a time depen- 
dent dielectric response e(t), the electromagnetic waves, 
V(x,t), are described by the wave equation: 



d 2 V(x,t) s d 2 e(t)V(x,t) 



dx 2 



dt 2 



(1) 



where So is the vacuum dielectric constant and c the 
speed of light in vacuum. By simple separation of vari- 



ables V(x,t) = ^2 k Vk(x)uk(t), the solution of this equa- 
tion may be expressed in the form: 

V(x,t) = J2<t)Mty kx = Y,Mt)e lkx - (2) 

k k 

After introducing the dimensionless variables, x — > x/L 
and t — > t/r, the equation for Uk(t) becomes: 



d 2 U k (t) 



dt 2 



k 2 s 2 f(t)U k (t) = 0, 



(3) 



where s = T cr / — 



and fit) = £ 1 (t). With e(t) periodic so 
that fit) is periodic, the solutions will be of the Floquet- 
Bloch jl3 type 



U k (t)=ex V (-iujt)U£(t), 



(4) 



which is the central property for photonic band-gap the- 
ory. 

By casting the problem in this form, it is reduced to de- 
termining the dispersion relation k = k(u). As is the case 
for traditional photonic crystals, the dispersion relation 
will contain the frequency regime where electromagnetic 
wave propagation is possible, as well as the "photonic 
gaps" where the propagation of electromagnetic waves is 
damped (or amplified) 13]. A possible realization of tem- 
poral dielectric modulation is to apply standing waves 
with time dependent intensity, in the form of one, two or 
three perpendicular laser beams, on a dielectric slab, as 
illustrated in Fig. [l] This concept is widely applied to 
create optical lattices within Bose-Einstein condensates 
[3 . Applying the beams on a spatially homogeneous or 
multi-dimensionally modulated dielectric slabs, a variety 
of space-time dielectric structures can be generated. 

First, we study the simple case of periodic dielectric 
modulation of a linear and homogeneous medium by 
modeling it (see Fig. (2J as a series of pulses character- 
ized by the two times r and ^. The parameter ^ is 
the duration of a single pulse of the dielectric function 
of the material e(t), and r is the period of repetition 
of the pulses. If the standing wave laser beam consists 
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where we have used the approximation 

dte-^-^U k (t)f(t)^U k (n) [ dte-^- n \f{t), 



which amounts to assuming the wave function Uk(t) to be 
constant during the pulse duration I/7 (see Fig. For 
the discrete set of amplitudes Uk(oo), when we integrate 
both sides of Eq. E3 and employ the relationship 



Ui 



duo 



: (m) = U k (t = m)= J ^-e*> 



FIG. 1: Illustration of concept and structure 




FIG. 2: Illustration of the modulation of the dielectric re- 
sponse 



we obtain the matrix equation: 

N 



U k {m) = ^2 ^n y mU kn (m), 



=1 



£ tempi 

2tt 



iuj{n—m) , .2 



(uo 2 - k 2 s 2 )(uo 2 + 7 2 ) 



(7) 



duo. (8) 



Without further approximations we find after some ma- 
nipulation (see supplement for details) that the disper- 
sion relation of this system can be expressed as: 



cos uo = h 2 temp - yj (hl emp ) 2 - 4,h° emp , (9) 



where 



h 1 = - 

,v temp 2 



cosh 7 + cos ks + rjtemp (cosh 7 + — sin ks) 

7 



of ultrashort pulses, a single pulse duration e(t) in the 
range 7 = sec can be realized. Similarly, a 

pulse separation cr comparable to the wavelength of the 
light in the medium is realistic. Specifically, we study 
the propagating waves in a 1 + 1 space-time scalar wave 
equation JJJ by assuming the dielectric constant to be of 
the form: 

N 

e(t) = 1+Y]g(t- mr) 

m=l 
N 

= 1+ J2 e -^flex P (- 7 \t-mr\). (5) 

m=l 

The specific choice (see Fig. El ) of the function g(t) 
is motivated by its convenient form for analytical treat- 
ment: For narrow pulse widths (7 — > 00), g(t) reduces to 
a f5-function and the model is then similar to the Kronig- 
Penney model. Applying this expression in Eq. © for 
the Fourier amplitude Uk(uo): 

- fcV) u k ( u ) = £ u k (m), (6) 

m—1 ^ ' ' 



h te = cosh 7 cos ks+r]temp (sinh 7 cos ks-\ cosh 7 sin ks) , 

7 

and 

7 3 

Vtemp- 2(^ + k 2 8 2 ) 6temP ' 

In Fig. 01 we illustrate this dispersion relation, Eq. 
The width of the gap is very small because of the realis- 
tic ^3 smallness of the modulation £te ™ p = 0.01 of the 
dielectric constant. One can easily see that the widths of 
the gaps are proportional to e terav . Figure 21 shows the 
width of the band-gap as a function of r and 7. This 
width is much larger for high frequencies but is relatively 
insensitive to 7. At this point it is instructive to note 
that for a traditional one dimensional photonic crystal, 
similarly defined by the dielectric constant 

e(x) = 1 + ^2 espat ^ Lx exp (-p\x - mL x \) , (10) 

m—l 

we obtain the dispersion relation by simply interchanging 
k with uo and ks with -. 
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4xl0 -1 ^ sec; 7 = 4; £ = 11; £t 
displayed in the upper inset. 



3: Spectrum of a "time photonic crystal" with r 
12 _ 11 o Q1 The "band-gap" 



Turning now to the more general question of multidi- 
mensional periodicity, we assume the dielectric constant 
e(r,t) to be a periodic function of both space and time, 
and solve the Maxwell equations for the electric field 
E(r,t). Again we have the Bloch type solutions with 
respect to space and time variables and so we introduce 
Bloch-Floquet parameters. The final dispersion relation 
will then relate k and uo. This means that in the case 
of combined spatial and temporal periodicity the electro- 
magnetic excitations spectrum is described by separated 
points represented by (k,(j), all belonging to an equiva- 
lent four-dimensional cube denning the four-dimensional 
Brillouin zone. For simplicity, we consider the dielec- 
tric constant e(r,t) to be of the type e(r,t) = e(r)e(t), 
i.e. we assume independent spatial and temporal mod- 
ulations of the dielectric media. The general solution of 
the wave equation is then 



E(r,t) = J2Vn(r)U n (t), 



(11) 



where for V n (r) and U n (t) we have 

d 2 
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FIG. 4: Width of the band-gaps as a percent of the mid- 
gap frequency for three frequencies r — 10 -13 sec - the lower 
= 10 -14 sec - the middle position curve, and r = 
sec - the upper curve. Triangles correspond to 7 = 
5, and squares to 7 = 3. The inset presents band-gaps for 
microwave region. An accessible experimental setup: r = 
10 _8 ,sec 7 — 10, Stemp — 0.2£o,£o — 300 was simulated. 



curve, r 
5xl0- 15 



and 



with 



V x V x V n (r) = -^ 2 £(r)V n (r) 



(13) 



dt 2 



e(t)U n (t) = -u 2 s 2 U n (t) 



(12) 



V-e(r)V n (r) =0 

and s = r CT / — = t 1 - and v being constants. In order 
to simplify the problem, we consider the wave equation 
for the scalar wave amplitude V(x,t). In this case of 
periodicity in one space dimension and in time, we have 
to take into account that the separation constant — v 2 is 
the same in both equations JT2J) and ([T3|). 

We can now establish the dispersion relation between 
the normalized frequency uo and the normalized wave- 
number k. For simplicity we first examine the case of 
s = 1, where the temporal and spatial periodicities are 
equal, L t = L x . We find that the dispersion of such a 
media is dominated by a vacuum -like dispersion rela- 
tion i.e. uj sk. The spectrum in this case is given 
mostly by the diagonal of the unit square, resembling in 
this sense a homogeneous medium. The actual disper- 
sion depicted in Fig. [5] shows this overall behavior but, 
as can be clearly seen in the inset, the underlying struc- 
ture is very complex. In fact the dispersion relation is a 
dense set of discrete points clustering around two sepa- 
rate branches. From this case we can easily construct the 
overall structure of the dispersion relation resulting from 
rational values of s ^ 1. For an integer values of s > 1 
the relation uo ~ sk results in values of uo > 1 and so 
this part of the dispersive branch is folded back into the 
standard zone. Similarly, for integer values of s _1 > 1 
the roles of uo and k are reversed and the k values have 
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FIG. 5: Dispersion relation of the temporally and spatially 
modulated medium with with r = 3xl0 -13 sec; 7 = 3; 
£o n ; L x 3xl0" 5 m ; p = 3 ;s x = 0.01; e t 0.01 
s=l. The inset shows a blow-up illustrating the underlying 
complex structure of the dispersion relation. 



to be folded back into the standard zone. Therefore, for 
rational values of s = M/N, where M and N are in- 
tegers leads to the repetition of M branches along the 
k axis and N branches along the uo axis. This is illus- 
trated in Fig. where the 5 = 2/3 case is displayed. 
When s is not a rational number, every single line de- 
termining the dispersion relations has a definite width. 
This means that in these cases any uo corresponds to an 
infinite number of k values, and vice-versa. This is com- 
pletely new situation, which has no analog in photonic or 
electronic band structures. The reason for this broaden- 
ing is the existence of small but finite band-gaps. The 
folding of the dispersion curves leads to a slight displace- 
ment of the equivalent curves. The largest displacements 
appear for the intervals of v where gaps are large. Then, 
at the larger values of z/, the curves become closer and 
when v tends to infinity they approach the correspond- 
ing straight line branches. This means that, when s is 
not a rational number, incommensurability between the 
spatial and temporal periodicities exists. In this case the 
dispersion relation consists of separated points and it is 
intriguingly complex, as is illustrated for s = y/2 in the 
inset of Fig. El Finally, we emphasize that all these ef- 
fects exist when the time variation of e(t) is rapid. When 
e(t) is a slowly varying function of time compared to the 
amplitude Uk(t), we can simply ignore its time depen- 
dence m. 

In summary, we have studied the effects of time vari- 
ation of the dielectric constant of different electromag- 
netic media, which leads to the existence of band-gap 
like phenomena. The physical reason for these effects 
is the necessity of synchronization of the phases of the 
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FIG. 6: Spectrum of the "space-time photonic crystal" with 
r = 3xl0~ 13 sec; 7 = 3; s = ll;s x = 0.01; e t = 0.01 
s = |. In the inset we represent a point spectrum - the 
incommensurability case of two basic lengths, where s = y/2 
and r = 3xl0 -13 sec. 

propagating waves with the external time periodicity of 
the medium. In cases of simultaneous space and time 
periodicity, we found an electromagnetic wave spectrum 
essentially described in an equivalent two-dimensional, 
three-dimensional or even four-dimensional cube denning 
the Brillouin zone. The folding of the dispersion relation 
when s = M/N ', with M and N integer, and broaden- 
ing of dispersion curves have been demonstrated. When 
s 7^ M/N the effects of incommensurability of the two 
internal lengths ( space L x and time L t = ^= ) exist, 
and point like spectra are exhibited. Modulation of the 
photonic band structure is suggested. 
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